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Abstract. - We present a generic formalism to describe Brownian motion of particles with 
intrinsic asymmetry and give predictions for the drift behavior in unbiased time-dependent force 
fields. Our findings are supported by molecular dynamics simulations. 



Introduction. — In view of their potential applica- 
tions in bio- and nanotechnology, Brownian motors have 
in recent years been the object of intensive research [1-4]. 
As is well known, the rectification of the thermal motion 
of Brownian particles involves the breaking of underlying 
symmetries. On the one hand, the system has to operate 
under nonequilibrium conditions to break the microscopic 
equilibrium symmetry of detailed balance. Spatial symme- 
1 try on the other hand is usually broken by applying asym- 
metric external forcing. The two most cited paradigms 
in this context are the flashing and rocking ratchets (see, 
e.g., [5]), in which an external space-periodic but space- 
asymmetric forcing using a ratchet-like potential is ap- 
plied. Somewhat surprisingly, the case in which an inher- 
ent asymmetry of the Brownian particle itself provides the 
spatial asymmetry has not been discussed in the context 
of periodic forcing. We will refer to such Brownian motors 
as intrinsic ratchets. In this letter, we will introduce and 
solve the equations of motion that generically describe this 
type of thermal rectification. 

i 

Generic equations of motion. The motion of a 
Brownian particle (speed v and mass M) is usually de- 
scribed by the following Langevin-Newton equation: 



dv 

M— = -jv + F + £. 
dt 



(1) 



Here 7 is the friction coefficient, F is an applied exter- 
nal force and £ a Gaussian white noise, whose intensity is 
determined by the fluctuation-dissipation relation. Equiv- 
alently, one can write the Fokker-Planck equation for the 



probability distribution P(v) of the speed, namely, 
7 F jk B T , 



d t P = d v 



M M 



M 2 



d v )P, 



(2) 



where ks is Boltzmann's constant and T is the temper- 
ature of the bath particles. The above equations can be 
derived from a microscopic description by assuming that 
the mass of the Brownian particle is much larger than that 
of surrounding particles. They form the starting point for 
deriving the properties of flashing or rocking ratchets. In 
fact, since one needs to apply spatially asymmetric forcing 
in these systems, the analysis is typically carried out at the 
simpler level of overdamped motion. The latter provides 
a closed description in terms of the position variable only 
and is known to be a very good approximation in most 
situations. As we will see below, we however do not need 
spatially dependent forcing for the rectification in intrin- 
sic ratchets. This greatly simplifies the analysis, even at 
the underdamped level. Indeed, when the forcing F is po- 
sition independent, the stochastic variable v is Gaussian, 
and it suffices to study the equations of motion for the 
first two moments of the velocity. By choosing as units 
of time, velocity and force the relaxation time t t — M/7, 
the thermal speed vt = yk^TJM and jvt, the follow- 
ing equations are obtained for the moments v\ = (v) and 



i'2 



1: 



dvi 

~dl 



= -vi + f, 



d «2 „ 



2M- 



(3) 
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We now argue that a minor modification of these equa- 
tions describes the case of intrinsic ratchets. We first note 
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that the possible asymmetry of the Brownian particle does 
not appear in the above equations, basically because the 
relaxation is described by linear response. As a result, the 
equation for the first moment, which is the central object 
of interest, is not coupled to the second moment. The 
asymmetry of the particle will appear at a next order of 
perturbation, at the level of nonlinear relaxation. Further- 
more, the resulting correction appearing in the equation 
for vi still has to vanish when operating under equilibrium 
conditions, i.e., when v 2 = 0. The simplest analytical cor- 
rection is thus a term of the form av 2 , where the constant 
a quantifies the strength of the asymmetry. Since this 
term acts like a perturbation on the first moment we can 
dismiss, to lowest order, the correction that will appear in 
the equation for v 2 ■ The intrinsic ratchet is thus described 
at lowest order (with a effectively playing the role of a 
small dimensionless parameter) by the following generic 
set of equations: 

-Vi + av 2 + f, 



dt 
dt 



In addition to the above handwaving arguments, we 
note that the equations of motion given in eq. (4) can 
be derived from microscopic theory of a Brownian particle 
moving in a bath of an ideal gas, by an expansion in the 
ratio of the mass of the gas particles (m) over the mass of 
the Brownian particle (M) [6-8]. Such a derivation also 
provides explicit expressions for the open parameters a 
and 7 (or r r ) behind eq. (4) in terms of microscopic quan- 
tities. Concrete examples for the cases of translational 
and rotational motion of an asymmetric object suspended 
in a thermalized gas will be given below. 

In the remainder of this letter, we focus on the recti- 
fication, i.e., the appearance of a non-zero average drift 
velocity, when the particle is subjected to an unbiased 
time-periodic force f(t). This scenario is the analogue of 
the rocking ratchet for particles with intrinsic asymmetry. 

Piecewise constant forcing. — Eq. (4) with time- 
periodic forcing f(t) has a mathematical structure similar 
to the Newton equation of motion for a parametric oscil- 
lator. It is therefore out of the question to find a general 
analytical solution. Instead we turn to the case of piece- 
wise constant forcing (square wave profile), viz., 

0<i<r/2: f(t)=f , 

t<T/2<r: f(t) = -f a . 

Introducing the vector notation 



(5) 



V(t) 



(*), 



(6) 



one readily finds the solutions in the separate time 
regimes, with f = fo and / = —fo respectively: 



v+(t) 

V-(t) 



A + (t)C+ 



(7) 
(8) 




Time 

Fig. 1: Time evolution of the first moment of the velocity v\ (in 
thermal speed vt units) subject to a modulating square force 
with periods 1, 4, and 16 (dashed, solid, dotted curve). The 
relaxation time r r of the particle is the time unit. The particle's 
asymmetry is a — 0.3 and force amplitude is fo = 0.4 (in units 
Mvt/t t , with M the particle's mass). 

The time-propagators A± are given by: 



A ± (t) 



j_ l-d± - 



4/o 



(3+d±)t/2 



-(3+d±)t/2 



4/o 



(3-d ± )t/2 



-(3-d±)t/2 



(9) 



(4) with d± = yT± 8a fo. B± are the steady state solutions: 



B± 



±/o/(l T af ) 
/oV(lT«/c) 



(10) 



The vector constants C± are specified by the assump- 
tion that we operate in the steady state regime, hence the 
velocity moments in eqs. (7) and (8) satisfy time-periodic 
boundary conditions, V+(0) = VL(r). This, together with 
continuity at t = t/2, V+(t/2) = V r _(r/2), leads to a 
solution, 

C+ = {A^ 1 A 1 - A- A l Az)-\Az l - A?)(B- - B + ), 

(11) 

C_ = (A^A 2 - A^A A )-\A^ 1 - A^)(B + - B_), 

(12) 

with At = A+(r/2),A 2 = A_(r/2),A 3 = A+(0),A 4 = 
A- (r) . We will not reproduce here the resulting expres- 
sion for the time-dependent average speed v±. It is ex- 
tremely cumbersome, and, strictly speaking, only valid to 
lowest order in the asymmetry contribution a. For an il- 
lustration of the typical time dependence of v\ we refer to 
fig. 1. We remark that our theory yields accurate results 
when afo <C 1. Under the described ratchet operation, 
the speed Vi and second moment v 2 are then of order fo 
and fo respectively. This means that the nonlinear cor- 
rection term av 2 in the equations of motion [eq. (4)] is a 
factor afo smaller than the linear v± term. 

The quantity of central interest is the resulting time- 
average net speed, being the average net displacement over 
a period r divided by this period: 



Vnct 



Vi(t) dt. 



(13) 



Again, the exact expression for v nct is extremely long. In 
any case, our approach is limited to small a, so it suffices 
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Fig. 2: Time-average net velocity v n ct as a function of (a) 
the asymmetry a of the Brownian particle, (b) the amplitude 
/o of the modulating force, and (c) the modulation period 
t. Exact analytical solutions for square modulation - solid 
curves - coincide, for the presented parameter range, with the 
first order approximation [eq. (14)]. Dashed curves represent 
results for harmonic modulation (upper curves) and sawtooth 
modulation (lower curves), eqs. (22) and (24), obtained from 
first order perturbation theory. Numerical integration results 
can not be distinguished from the analytical solutions in the 
shown graphs. If not in the abscissa, parameter values are 
q = 0.3, /o = 0.4, t — 1. Units for velocity, time and force are 
the particle's thermal speed vt, relaxation time r r and Mvt /t t 
(mass M) respectively; a is dimensionless. 



to reproduce the lowest order term in a: 



4 T 

Vnct ^ a/o (1 tanh - 

T 4 



(14) 



We note that the next term in the expansion in a is an 
order of magnitude (a/o) 2 smaller. 

We make the following observations. First, there is 
no directed motion, v nct = 0, in the absence of forcing, 
/o = 0, or when the particle has intrinsic symmetry, a = 0. 
Second, v nc t is an uneven function of a, hence an inver- 
sion of the asymmetry, a — ► —a, results in an inversion of 
the speed of net motion. We represent v nct as a function 
of the asymmetry, a, the amplitude of the force, fo, and 



the period, r, in fig. 2 (solid curves). In all three cases, 
the lowest order approximation, eq. (14), is in fact indis- 
tinguishable from the exact result for the chosen range of 
values of a and Jo- Finally we note that the maximum 
speed ujj™ ~ a/p is reached in the limit of very slow mod- 
ulation, r — > oo. Since this speed is expressed in units of 
thermal velocity, we conclude that one can reach high net 
drift speeds, comparable to thermal speeds, by applying 
unbiased periodic forcing of small to moderate intensity 
to intrinsic ratchets. 

Other periodic forcings. — To investigate the ef- 
fect of other types of periodic forcing, such as harmonic or 
symmetric sawtooth [cf. fig. 3(b)], we resort to a perturba- 
tional solution of eq. (4). As the contribution due to the 
intrinsic asymmetry a is considered small, we can make 
the following first-order ansatz for the velocity moments: 



V2 = "2,0 + OiV 2 ,l ■ 



(15) 
(16) 



With this ansatz, eq. (4) can be solved to first order in a 
for arbitrary periodic force fields f(t), yielding the steady 
state solutions 



w i,o0) 

«l,l(*) 
W2,l(*) 



— OO 

t 



— oo 

t 



— oo 

t 



di'e-t*" 4 ') /(*'), (17) 

dt ' e -2(t-t') 2 /(i> 1>0 (t'), (18) 

di'e-C-'VoCi'), (19) 

di , e -2( t -0 2 /(t') Vll (^). (20) 



Using f(t + t) = f(t) it is easy to show that these ex- 
pressions are indeed periodic with periodicity r. The re- 
sults for the time evolution of the first moment v% under 
harmonic or sawtooth forcing are reproduced in fig. 3(a); 
they are indistinguishable from numerically integrated so- 
lutions of the original eq. (4). For comparison, the results 
for the square wave profile are also included. 

With regard to the net velocity u nct as defined in 
eq. (13), we observe that for unbiased symmetric forc- 
ings /(i + r/2) = -/(*), and thus /J" At' v lfi {t') = [see 
eq. (17)], so that v nct is given by 



V net = CUT 



d*Vi(0 



(21) 



Again, w net is zero when a = 0, consistent with the notion 
that no directed motion occurs for symmetrical particles. 

It is straightforward to recover the net velocity for 
square forcing to lowest order in a, eq. (14), from (21). 
Similarly, for harmonic driving, f(t) = f sin(27rf/r), we 
find a time-average net velocity 



"not 



«/o 2 



2 4tt 2 + t 2 



(22) 
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Fig. 3: Time evolution of (a) the first moment of the veloc- 
ity vi (in thermal speed vt units) for the periodic forcings in 
(b) (correspondence by line style). For harmonic and symmet- 
ric sawtooth profiles the first order approximation [eq. (15)] 
is used, for square modulation the full analytical solution is 
shown. In all three numerical solution of the dynamic 

equation is indistinguishable from the analytical results in the 
graph. Parameter values are period r = 4, force amplitude 
/o = 0.4 and asymmetry a = 0.3. Unit of time is the parti- 
cle's relaxation time r r , unit of force is Mvt/t t , with M the 
particle's mass. 



For sawtooth forcing, 

0<t<r/2: f{t) = f (t - r/4)/(r/4) , 
t< T /2<T: f(t) = /o (3r/4 - t)/( T /4) . 



a net speed 



Vnct 



<xfo 
3 



:-! | ij tanh^ 



(23) 



(24) 



is obtained. These first order results for the net velocity 
v nc t are compared with the analytical solution for square 
forcing in fig. 2 (dashed curves). We conclude that the 
resulting drift behavior is very similar in all three cases. 
In fact, comparing eqs. (14), (22) and (24), we see that the 
differences become very small, and even vanish for slow 
forcing t — * oo, if, instead of using the same amplitude 
for the three modulations, one considers the same average 
quadratic amplitude, i.e., if one replaces /o/v2 — ► /o in 
eq. (22) and f /V3 -> fo in eq. (24). 

Microscopic models. — As already mentioned, the 
structure of eq. (4) can be obtained from kinetic theory 
of microscopic models, that describe a small, non-trivially 
shaped object (mass M) moving under the influence of 
collisions with a surrounding bath of gas particles (mass 
m), by an expansion in the mass ratio e = y/m/M [6-8]. 
Such a procedure also provides explicit expressions for the 
parameters r r and a of the intrinsic ratchet. 

For a three-dimensional asymmetric object of arbitrary 
convex shape, that is confined to move along a fixed z-axis 
(translational motion), one obtains [8] 



a — \J 7r/8e 3 r r cr3, 



(25) 
(26) 





Fig. 4: Idealized realizations of the intrinsic ratchet, (a) For 
translational motion: a conical shape with axis along the free 
direction of motion. Indicated is the resulting sense of net 
motion under unbiased forcing, (b) For rotational motion: a 
right triangular prism, with three suggested locations for the 
rotation axis. The resulting net rotation sense under unbiased 
forcing is indicated. 



where the geometry dependent moments o~ n are given by 



™ Js 



(27) 



with p being the particle density of the gas, and where the 
integral is over the surface of the asymmetric object. e±\ z 
is the component in the free direction of motion (z) of the 
outward unit normal vector e± at its surface. 

To get an idea of the actual net velocity a realistic setup 
of an intrinsic ratchet can attain, we use eqs. (14), (25)- 
(27) to calculate the speed for a cone-shaped silica (Si02) 
Brownian particle with half opening angle 30° and 10 nm 
base diameter [see fig. 4(a) for a schematic representation]. 
The cone's axis is along the free direction of motion. In 
air, the ratio M/m is about 10 000 and the asymmetry 
parameter is a = 0.003. At temperature T = 300 K the 
relaxation time is r r = 7.5 ns. An amplitude /o = 10 of 
unbiased square forcing then corresponds to 1.9 pN in real 
units and is well within the accuracy range of our theory: 
ctfo = 0.03. These conditions produce a maximum speed 
of 

^net = 0.88 m/s or 30% of the thermal speed of the 
particle. Note that the direction of the particle's motion 
is towards the apex of the cone. 

Rotational Brownian motion. For simplicity of 
presentation, we started with the generic equations of mo- 
tion for one-dimensional translational Brownian motion 
of an asymmetric particle. In practice, this supposes that 
the particle is constrained to move on a track. The so-far 
presented discussion of the intrinsic ratchet can however 
be repeated, with minor modifications, for the rotational 
Brownian motion of chiral objects, with angular veloc- 
ity u! and moment of inertia /. With an adaptation of 
the expressions for the relaxation time r r =7/7 and the 
thermal velocity vt = \JksT /I as units of time and an- 
gular velocity, and / now signifying a torque, this leads 
to the same generic equations of motion [eq. (4)] for the 
moments v\ — (w) and v-i — (u 2 ) — 1. Microscopic the- 
ory [9, 10] yields the same expressions for the relaxation 
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time 7> [eq. (25)] and asymmetry coefficient [cq. (26)], but 
now with geometrical moments given by 



<?n = pJ^L I dS [(ej_ x e r )Uf , (28) 
V 7rm J s 

where the axis of rotation is taken to be parallel to the 
z-axis. Again, the integral is over the surface of the ob- 
ject and e_L is the outward unit normal vector on the sur- 
face. e r is given by r/rg, with r denoting the position of 
a surface element measured from the axis of rotation (the 
z-component in r is irrelevant), and r = \J I /M being 
the radius of gyration of the object. Due to the chosen 
orientation of the rotation axis, only the z-component of 
e± x e r appears in the expression for a n . With these new 
notations and units, the previous results, in particular the 
expressions for the time-average net velocity [eqs. (14), 
(22), (24)], remain valid. 

In view of the technological potential of the rotational 
setup, and in order to get an idea of the order of magni- 
tudes involved, we again consider a realistic physical real- 
ization. A silica triangular prism of height 10 nm and with 
right triangular top and bottom surfaces (sides: 10 nm, 
10 nm, 14 nm) is connected with a rotation axis at one 
of its vertical edges, cf. fig. 4(b). Operating in air, the 
ratio M/m is about 23 000. If the axis is connected to 
the 90° corner edge, a = and no rectification or net 
rotation will occur. Connected to the 45° corner edges, 
a = ±0.0016 and relaxation time is r r = 8.1ns, at air 
temperature T — 300 K. A torque amplitude (for square 
forcing) of /o = 10 corresponds to 2.2 x 10 _20 Nm and 
produces a maximum net rotation frequency of 6 MHz, 
16% of the thermal frequency. 

Molecular dynamics simulations. In the follow- 
ing we verify our generic theory for the intrinsic ratchet 
with molecular dynamics simulations. As concrete micro- 
scopic system, we consider the prism from fig. 4(b) sur- 
rounded by a thermalized bath of ideal gas particles, and 
perform event-driven simulations of its rotational Brow- 
nian motion. The rotation axis is chosen to be located 
slightly [4nm in the units of fig. 4(b)] outside the prism 
in the plane given by one of the prism surfaces merging at 
the 90° edge, and is oriented parallel to this edge. 

Exploiting the homogeneity of the prism along the di- 
rection of the rotation axis, the simulations are carried out 
in a (projected) two-dimensional space, where the ratchet 
object is given by the right triangular top (or bottom) sur- 
face of the prism, and where the rotation axis is reduced to 
a point-like center of rotation. The rotation center is posi- 
tioned at the center of a quadratic box, containing an ideal 
gas of point particles (mass m). The box walls 'absorb' 
gas particles upon collision, but also randomly 'emit' new 
particles (into the box' interior) such that the gas prop- 
erties, in particular density p and Maxwellian equilibrium 
distribution, are preserved. In this way, an infinitely large 
reservoir of gas particles is realized. 
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Fig. 5: Molecular dynamics simulation results (dots with error 
bars) compared with theory (curves) for the rotational intrinsic 
ratchet. Time-average net velocity v nct is shown as a function 
of the modulation period r of zero-average square forcing. The 
simulation results are obtained from averaging over typically 
5000 realizations (20 000 for the smallest driving periods) with 
about 17 000 collisions performed per realization, correspond- 
ing to a simulation time of 1000 periods for the fastest and 20 
periods for the slowest modulation; the error bars characterize 
the remaining statistical uncertainty. Simulation parameters 
are p = 0.25, M = 50, m = 1, side lengths (1, 1, y/2) of the 
right triangle, and amplitude 0.15 of the external square torque 
(see also main text), resulting in a — 0.0233 and fo = 1.32. 
Units are relaxation time r r (time), thermal angular velocity 
vt = y/ksT/I (angular velocity) and Ivr/r r (torque), with / 
the inertial moment. 

Collisions between gas particles and the triangle are de- 
tected by numerically solving the exact equations of mo- 
tion for the point in time of the impact. At each collision, 
the speed of the gas particle and the rotational velocity of 
the triangle are changed according to the rules for elastic 
collisions, neglecting tangential forces [9]. In between col- 
lisions the triangle is accelerated by an external constant 
or periodically switching torque (square profile). 

In fig. 5 the net rotational speed of the triangle under 
the action of a periodically switching torque is shown for 
different values of the modulation period r. The agree- 
ment between simulation results for the time-average net 
velocity w nct and theory, eq. (14), is excellent. We also 
compared the asymptotic rotation of the triangle under 
constant but opposite torques ('infinite' driving period r) 
with the theoretical result, eq. (10), and again found ex- 
cellent agreement. 

Stop-and-go motor. — We finally discuss an alter- 
native approach to intrinsic ratchets, anticipated in [11] 
and further worked out in more detail in [12]. The so- 
called stop-and-go motor consists of an asymmetric parti- 
cle which is periodically stopped, for example by an array 
of traps or binding sites that can be activated or deacti- 
vated at will. The basic assumption is that the (thermal) 
energy of the Brownian particle is changed when it is sub- 
jected to the trapping mechanism. This energy exchange 
results in a specific value of the second velocity moment at 
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holds, even for large i>2, s - 
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-0.004 



12 3 4 

Stopping interval 

Fig. 6: Main figure: average (angular) net velocity v nc t of the 
rotating stop-and-go ratchet as a function of the time interval 
r s between stopping events. In the trapping mechanism the 
second moment is set to «2,s = 1 (upper curve), V2,s = — 1 
(lower curve), with V2 = {w 2 } — 1. Insets: maximal average 
net velocity v nct (for optimal t s ) as a function of the second 
moment of the velocity «2, s in the trap. Lines correspond to 
theory, dots represent results from molecular dynamics simula- 
tions, using the same setup as for fig. 5 (but without external 
torque). The net speed is obtained from simulating 1000 in- 
dependent stopping intervals t„ per realization and averaging 
over 10 000 realizations for the three smallest stopping intervals 
and over 5000 realizations otherwise. Shown angular velocity 
is in units of thermal speed vt = y/ksT/I and time in units of 
relaxation time TV- The asymmetry parameter is a — 0.0233. 



the stopping sites, V2, s - For v 2 . s ^ the energy exchange 
in the traps induces a deviation from thermal equilibrium 
conditions, and this process will result in sustained di- 
rected motion with an average net speed v ne t, being the 
average distance traveled by the particle in a time interval 
t s between the stopping events, divided by r s . A simple 
analytical calculation, starting from eq. (4) with / = 0, 
gives an exact expression: 



Vnet 



aV2,s 
277, 



-t.\2 



(29) 



The sense of motion is determined by the sign of a and 
of V2, s , which is negative (positive) when the particle's 
thermal motion is reduced (enhanced) by the trapping 
mechanism. A stopping interval r° ~ 1.26, given by the 
solution of e Ts — 2r s + 1, yields a maximum net veloc- 
ity (wJJet* — 0.204 av2, s ) and an optimal distance between 
binding sites (~ 0.256 av2, s , expressed in units VTT r ). 

These theoretical predictions are confirmed in a molecu- 
lar dynamics simulation of the stop-and-go mechanism ap- 
plied to a rotating (chiral) object, using the setup based on 
the prism in fig. 4(b), as in the previous section. In fig. 6 
the resulting average net velocity v no t as a function of dif- 
ferent stop intervals r s , for the values V2, s = — 1 and 1 is 
shown. In the insets of fig. 6 we include the molecular dy- 
namics results for w net at the optimal stopping interval as 



Conclusion. — Intrinsic ratchets are characterized by 
an inherent asymmetry of the Brownian particle itself 
breaking the spatial symmetry. A generic formalism for 
the dynamical behavior enables us to quantify the net 
particle velocity under unbiased periodic forcing. Molec- 
ular dynamics simulations of a rotational setup confirm 
the validity of this formalism. We predict drift speeds 
comparable to thermal speeds for nanosized asymmetric 
Brownian particles under ratchet operation. The relative 
simplicity of the setup (one heat bath, external symmet- 
ric forcing) could open avenues to experimentally test the 
nonlinear contribution of intrinsic asymmetry crucial in 
this and other work [6-15]. 
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